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According to the recent studies, the gravitational wave (GW) echoes are expected to be generated by quark 
stars composed of ultrastiff quark matter. The ultrastiff equations of state (EOS) for quark matter were usually 
obtained either by a simple bag model with artificially assigned sound velocity or by employing interacting 
strange quark matter (SQM) depicted by simple reparameterization and rescaling. In this study, we investigate 
GW echoes with EOSs for SQM in the framework of the equivparticle model with density-dependent quark 
masses and pairing effects. We conclude that strange quark stars (SQSs) can be sufficiently compact to possess 
a photon sphere capable of generating GW echoes with frequencies in the range of approximately 20 kHz. 
However, SQSs cannot account for the observed 72 Hz signal in GW170817 event. Furthermore, we determined 
that quark-pairing effects play a crucial role in enabling SQSs to satisfy the necessary conditions for producing 


these types of echoes. 
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I. INTRODUCTION 


The detection of gravitational waves (GWs) by the LIGO 
and Virgo collaborations [1-5] has significant implications in 
astrophysics, opening up new avenues for research on black 
holes, neutron stars, and other dense stellar objects [6, 7]. It 
was suggested that GW echoes [8—10] were expected as a 
generic feature of quantum corrections at the horizon scale in 
the postmerger GW signals from binary coalescence events, 
particularly those involving black holes. In general, the emis- 
sion of GW echoes requires dense stellar objects featuring a 
photon sphere located at Rp = 3M, where M denotes the 
mass of the dense object [11-14]. The photon sphere can par- 
tially trap the GW to produce the GW echoes. 

This type of a photon sphere can be found both in black 
holes [15] and superdense stars [16, 17], whose radii should 
be smaller than R,. For black holes, besides the photon 
sphere, the GW echoes necessitate a secondary reflecting sur- 
face to circumvent GW absorption, which may be related to 
quantum effects in proximity to the black hole horizon [18]. 
The possible occurrence of GW echoes at a frequency of 
approximately 72 Hz, with a statistical significance level of 
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4.20, was first investigated in Ref. [19] for the GW170817 
event, where the GW echoes were interpreted as a result of 
quantum effects close to the black hole horizon. 


However, Ref. [20] indicated that GW echoes can be gen- 
erated not only by quantum corrections at the horizon scale 
but also by exotic super-compact objects [21, 22]. Moreover, 
a simplified incompressible equation of state was utilized, re- 
vealing that a highly compact stellar object with radius close 
to the Buchdahl’s radius [23] Rp = 9/4M is required to pro- 
duce GW echoes with such low frequencies. Specifically, the 
Buchdahl’s radius is the minimum radius for compact stars 
[23]. Additionally, in order for compact stars to produce GW 
echoes, the radii of compact stars must fall within the range 
of Rg < R < Rp. This compactness criterion essentially 
precludes the realistic equation of states (EOSs) for neutron 
stars [20]. Consequently, researchers are exploring alternative 
compact stellar objects, such as strange quark stars (SQSs) 
[24], which consist of strange quark matter (SQM) [25-28]. 


In Ref. [29], the authors verified whether a more realis- 
tic EOS of quark stars can emit GW echoes. They employed 
the confined-isospin-density-dependent-mass model with ad- 
ditional scalar and vector Coulomb terms of SQM [30] and 
confirmed that SQSs with realistic EOS cannot be catego- 
rized as ultra-compact objects to feature a photon sphere to 
generate GW echoes. Moreover, in Ref. [31], the authors 
utilized an interacting quark matter EOS unifying interacting 
phases via a simple reparameterization and rescaling. They 
found that GW echoes are possible for QSs with large cen- 
ter pressure. Furthermore, GW echoes were examined in 
f(R,T) gravity metric formalism within MIT bag model and 
the color-flavor-locked (CFL) EOSs. The authors indicated 
that, under some considerations, the realistic interacting QM 
can lead to stellar structures, which are sufficiently compact 
to feature a photon sphere outside the stellar boundary, and 
thereby, can echo GWs [32, 33]. Additionally, the author 


investigated the GW echoes from SQSs for various EOSs, 
including MIT bag model and linear and polytropic EOSs. 
However, only the MIT bag model and linear polytropic EOSs 
were found to emit GW echoes at a frequency range of ap- 
proximately tens of kilohertz [34]. In a recent study [35], 
the GW echoes produced by strangeon stars composed of 
strange-cluster matter, which is in the solid state, were in- 
vestigated. The authors recasted the EOS of strange-cluster 
matter into dimensionless forms via reparameterization and 
rescaling. Furthermore, they concluded that strangeon stars 
are typically compact enough to have a photon sphere. This 
sphere reflects the GWs that fall within the gravitational po- 
tential barrier, producing GW echoes with a minimum echo 
frequency of approxiamtely 8 kHz, extending even to fre- 
quencies as low as O(100) Hertz. 


Nonetheless, in the aforementioned studies, the investiga- 
tions concerning GW echoes have primarily concentrated on 
either a basic MIT bag model or parameterized EOSs for 
SQM. In this study, we employ EOSs for SQM based on the 
equivparticle model with density-dependent quark masses. 
Within this model, quark masses are scaled according to the 
baryon number density, replicating the complex interactions 
among quarks [36]. 


This model was initially developed as a quark mass- 
density-dependent model [37—40], and later renamed as the 
equivparticle model after explicitly introducing the concept 
of effective quark chemical potential [41]. Given its clear 
physical picture and accurate thermodynamic treatments, the 
model has been widely utilized in the study of quark matter 
properties and structures of quark stars [42—46]. In Ref. [47], 
we investigated the symmetry energy of SQM and tidal de- 
formability of SQSs within this model. Our findings indicate 
that the region of absolute stability for SQM can be signifi- 
cantly widened with sufficiently large isospin symmetric pa- 
rameter, yielding results that simultaneously satisfy the con- 
straints imposed by astrophysical observations of PSR J1614- 
2230, with 1.928 + 0.017 Mọ [48] and tidal deformability 
70 < A4 < 580 measured in the GW170817 event [7]. 
On application of the EOSs in this model, we determine that 
SQSs can produce GW echoes with the frequencies in the or- 
der of kilohertz, which is consistent with previous findings, 
provided that they are composed of SQM in the CFL phase 
[32, 33]. 


The paper is organized as follows: Section II provides a 
brief overview of the EOSs for SQM in CFL phase within 
the equivparticle model. Section III presents the process of 
calculating the GW echoes and scrutinizes the corresponding 
numerical results. Finally, a concise summary is presented in 
Sect. IV. 


II. EOS OF CFL QUARK MATTER IN EQUIVPARTICLE 
MODEL 


The thermodynamic potential density for CFL SQM can be 
expressed as: 


Qcorn = = E f(m- p yee end) 9? dp 
T 
i= 2 
3A? p? 
= Ep, (1) 
where gi = 2 x 3 = 6 corresponds to degenerate 
factor, v denotes the common Fermi momentum, ñ = 


(Hu + Ha + Hs) /3 is the averaged chemical potential of 
quarks, A denotes the quark pair energy gap, B denotes the 
bag constant, which can be expressed as B!/4 = 180 MeV in 
the following calculation, and m; denotes density-dependent 
quark mass which can be scaled as: 
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where my) = Mao = 5 MeV and m.o = 100 MeV are the 
quark current masses [49], my denotes the interacting part 
and is the same for all flavors, and D denotes a model pa- 
rameter representing the strength of confinement. Due to the 
lack of accurate value of A, in this model A is considered 
to be a free model parameter. It should be emphasized that 
the quark chemical potential u; is effective. Furthermore, the 
real chemical potential 4; rea can be related to the effective 
one due to the density-dependent quark mass [41]. To ensure 
maximum paring, the common Fermi momentum v can be 
obtained by taking the derivative of Ocrr with respect to v, 
i.e., OAcr,/Ov = 0, yielding 


So 4m? = 3). (3) 


The integration in Eq. (1) can be conducted as follows: 
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from which one can readily obtain the quark number den- 
sity for each quark flavor according to the relation n; = 
—0NcrL /ðpi, i.e., 
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(v? + 2A7p)/x?, (5) 


where np = (Nu t+na+ns)/3 denotes baryon number density. 
The energy density is as follows: 


E = Qcort + 5 Lini = NcFL + 3ni (6) 


Due to the density-dependent quark mass, the pressure of 
the system is: 


OQcrL dm; 
Om, dn’ 
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where the second term in the right side of the equal sign is 
crucial to guarantee the thermodynamic consistency, and the 
derivatives are, respectively, as follows: 
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HI. GW ECHO FREQUENCY OF CFL SQS IN 
EQUIVPARTICLE MODEL 


To generate the GW echoes, the EOSs of SQM should be 
sufficiently stiff to feature a photon sphere at Rp = 3GM. To 
examine how the stiffness of EOS changes with model param- 
eter D and A, we express the velocity of sound as function of 
VD (upper panel) and A (lower panel) with respect to differ- 
ent baryon number densities provided in Fig. 1. Based on the 
upper panel in Fig. 1, at a fixed A = 150 MeV, the velocity 
of sound decreases as VD increases, indicating a softening 
of the EOS with larger D. Conversely, with fixed VD = 100 
MeV, the velocity of sound increases with A, which implies 
that the EOS becomes stiff for large A. Therefore, based on 
Fig. 1, EOS can be inferred to become more stiff by simulta- 
neously reducing D and increasing A. Hence, a stiff EOS can 
be realized by setting D = 0 and simultaneously adopting a 
large A. However, the value of A cannot be arbitrarily chosen 
for a given D and is dependent on the density of SQM. To in- 
vestigate this point, the following steps should be considered. 

First, by combing Eqs. (3) and (5) and then eliminating 
the averaged quark chemical potential ñ, >; \/v2 + m? = 

3. (7? — v?) can be obtained. Based on this, a function 


2X2 ; ; 
of v can be introduced, i.e., 


BY) = sep(n'm—9)— oF + mB, (v2 0). 1 


Here, note that the purpose of introducing the function ®(v) 
is only to determine the relation between model parameters 
D and A. Additionally, v should be considered as a variable 
of function ®(v), although it has the meaning of Fermi mo- 
mentum. 

Then, the derivative of ®(v) can be easily obtained with 
respect to v via d®(v) /dy, 
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Fig. 1. Velocity of sound as a function of model parameters D (up- 
per panel) and A (lower panel) with different given baryon number 
densities np. 


where ae) is observed to be consistently negative (equal to 


zero only when v = 0). This suggests that ®(v) increases 
with decreasing v, and thereby, indicating the existence of a 


maximum value of ®(v) at v = 0, namely 


37? np 
®max(v = 0) = IS a (12) 
For Eq. (12), we should indicate that v = 0 is only a limiting 
case, and its purpose is to provide the D-A window. Addi- 
tionally, in the following calculations, non-physical values of 
v, such as v = 0, will not be considered. 

If Eq. (10) admits a solution for v, then it implies that the 
maximum value ®,,.x( = 0) should be no less than 0. This 
in turn yields Eq. (12), i.e., ®max(v = 0) > 0, which leads 
to inequality 


3Np 
230, mi(m, D)’ 


where the quark mass m; is explicitly expressed as a function 
of baryon number density np and model parameter D. 

Evidently, by assigning a value to the model parameter D, 
the maximum value of A can be expressed as a function of 
baryon number density. Furthermore, the maximum value 
of A increases with n, for a fixed D. This implies that the 
highest attainable value of A is contingent upon the potential 
minimal value of the density of SQM for a specific D. Con- 
sidering that the saturation density of normal nuclear matter 
is no © 0.16 fm~ 3, we assume that the density of SQM is not 
smaller than no, and the minimum baryon number density of 
SQM is designated to be np min. To establish a range of val- 
ues for A, we set Np min = No, 2no, and 3no in Eq. (13). 

In Fig. 2, we present the model parameter window in the 
A—D'/? diagram for different minimum baryon number den- 
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Fig. 2. (Color online) Model parameter window A — D'/? with re- 
spect to the given lowest baryon number density np, min for SQM. 
The region bounded by the solid, dashed, and dotted black lines cor- 
respond to the solvable EOSs for SQM with the minimum densities 
that are not less than no, 2no, and 3no, respectively. The red line 
illustrates the correlation between A and D\/ 2, which can lead to 
the most massive SQSs, located precisely on the photon sphere line. 
The typical model parameters are denoted by solid squares and la- 
beled as A, B, and C. 


sities, nb min, Of SQM. For instance, if the model parameters 
lie within the region below the solid black line, a solvable 
EOS for SQM with a minimum density of Np min = No can be 
available. The black dashed and dotted lines serve the same 
purpose but correspond to minimum densities of 2n9 and 3noọ, 
respectively. Referring to previous study on CFL SQM within 
this model [39], when assigning a value of A = 100 MeV, we 
find that a solvable EOS spans a broad range of D1/?, from 
0 to beyond 130 MeV. Notably, as D rises, the maximum 
value of A decreases for a specified np min. Additionally, in- 
creasing Nb min Significantly expands the parameter window. 
For subsequent calculations, we select three representative 
model parameter sets: (V D/MeV, A/MeV) = (70,100), 
(70, 359.1), and (70,500), denoted as sets A, B, and C, re- 
spectively, marked by solid squares in Fig. 2. 


Before delving into the mass-radius relation for SQSs, we 
m 


should first check whether condition A > Z can be fulfilled. 
This guarantees that the CFL SQM can stably exist [50]. In 
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Fig. 3, we show the difference between A and rr 
baryon number density ranging from no to 10no for the typi- 


cal model parameters chosen in FIG. 2. Based on this figure, 
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TA is positive, which implies 


that the condition A > me is fulfilled. Additionally, with in- 
creasing baryon number density, the difference becomes sig- 
nificant. Furthermore, the large A provides large difference 


the difference between A and 
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To explore whether the SQSs in the current model can pos- 
sess a photon sphere and subsequently produce GW echoes, 


one should first calculate the EOSs for CFL SQM. Following 
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Fig. 3. Difference between A and Zu Versus baryon number den- 
sity. The difference is observed to be positive, which implies that the 


condition A > z is fulfilled. 
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Fig. 4. Mass-radius relation of SQSs in the EOSs within the equiv- 
particle model with typical model parameter sets labeled as A, B, 
and C. The values of parameter sets are (V/D/MeV, A/MeV) = 
(70, 100), (70, 359.1), and (70, 500) for A, B, and C, respectively, 
as shown in Fig. 2 with solid squares. The maximum masses in 
these cases are indicated with solid dots. 


this, the mass-radius relation for hydrostatically equilibrated 
SQSs can be determined by solving the TOV equations. The 
TOV equations are: 


dP G(m + 4rr? P) 

dr err) 2Gmr—r? ’ (14) 
Cs eg (15) 
dr 


where Æ and P denote the energy density and pressure at 
radius r, G = 6.707 x 10745 MeV~? is the gravitational 
constant, and m(r) denotes the gravitational mass within the 
radius r. 

The resulting mass-radius relations are illustrated in Fig. 
4. Based on this figure, when A = 100 MeV, the maximum 


mass for case A is determined to be below the photon sphere 
line. This implies that no GW echoes can be generated in this 
case. While when A = 359.1 MeV, the maximum mass of 
the SQS nicely locates at the photon sphere line, signifying 
the critical state for the SQS to possess a photon sphere. For 
a larger value of A = 500 MeV, the EOS corresponding to 
the maximum mass of SQSs surpasses the photon sphere line, 
indicating its capability to produce GW echoes. This suggests 
that the configurations for the SQSs, ranging from point 1 to 
point 2 in the detailed review profile shown in Fig. 4, ex- 
hibit the potential to generate GW echoes. Thus, for a fixed 
value of VD (specifically, VD = 70 MeV in this instance), 
the maximum mass of SQS, represented by solid dots in Fig. 
4, increases with A. This correlation is further emphasized 
in Fig. 1, where the velocity of sound rises concomitantly 
with A. Moreover, for a constant VD, a specific A value can 
be determined that positions the maximum mass of the SQS 
exactly on the photon sphere line. 

In Fig. 2, the red line illustrates the correlation between A 
and D!/2, which can lead to the most massive SQSs located 
precisely on the photon sphere line. Furthermore, with an 
error in the level of 0.1%, it can be fitted as 


3 a 
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where the coefficients are cj ~ 114.5837, co ~ 85.3983, 
and c3 ~ 52.5021. Based on this equation, parameter A 
is found to increases with VD and A, and the first term 
Amin © 355.6398 MeV is the minimum value of A at D = 0. 
This means that the SQM in the CFL state within the cur- 
rent model can satisfy the requirement of featuring a photon 
sphere for SQS. Consequently, the parameter sets (V D, A) 
that qualify SQSs to produce GW echoes should lie above the 
red line illustrated in Fig. 2. Hence, only the parameter set 
C aligns with the SQS configurations requisite for generat- 
ing GW echoes. Additionally, from Eqs. (13) and (16), we 
deduce that the least density necessary for SQM to feature a 
photon sphere within this model is around np min œ 0.1223 
fm-°. 

To derive the frequency of GW echoes, the time taken for 
light to traverse from the center of the star to its photon sphere 
[20] should be computed, i.e., 
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When 0 < r < R, m(r) = Z pr’, The gravitational po- 
tential, (r), can be derived by integrating the subsequent 
differential equation, 1.e., 


dð(r) 1 dP 


— 1 
dr E+P dr’ te) 


which can be solved together with the TOV equations in Eqs. 
(14) and (15); When R < r < 3GM, m(r) = M, e?® = 
1 — 2M/r, and P = p = 0. Here, P and E denote the 
pressure and energy density, respectively, as determined by 
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Fig. 5. Frequency of GW echo wecho and the mass of SQS as func- 
tions of central pressure Po. 
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Fig. 6. (Color online) GW echo frequency wecho for the most mas- 
sive SQS in the model parameter window. 


the EOS. Finally, the GW echo frequency can be evaluated 
by Wecho = T j Techo» 

In Fig. 5, we present the GW echo frequency wecho and the 
mass of SQS versus central pressure Po for case C from point 
1 to point 2 in the detailed review profile in Fig. 4. Evidently, 
as shown in Fig. 5, the mass of SQS increases with central 
pressure, while the GW echo frequency decreases. In fact, 
this can be interpreted as follows: with higher pressure, the 
SQS will have a larger mass. Based on Rp = 3G'M, the pho- 
ton sphere will be situated further from the center of the SQS. 
This will result in light requiring more time to travel from the 
center of the SQS to the photon sphere, thereby leading to 
a smaller Wecho. Therefore, given the appropriate model pa- 
rameters, the most massive SQS exhibits the lowest GW echo 
frequency. 

In Fig. 6, we show the GW echo frequencies for the most 
massive SQSs within the model parameter window where A 
and vD are in the range of 400-800 MeV and 0-130 MeV, 
respectively. According to Fig. 2, the selected model param- 
eter ranges ensure that the most massive SQSs are positioned 
above the photon sphere line, thereby capacitating them to 


generate GW echoes. Notably, with a constant A, a negli- 
gible shift is observed in the GW echo frequency as D es- 
calates. Conversely, with a stable D, the GW echo frequency 
undergoes notable variation with changes in A; specifically, it 
transitions from approximately 22.8 kHz at A = 400 MeV to 
approximately 18.8 kHz at A = 800 MeV. In conclusion, the 
GW echo frequencies predominantly oscillate around 20 kHz 
in this model, which aligns with the GW frequency magni- 
tudes reported in prior studies [20, 24, 31]. Additionally, the 
unmarked region designated as imbalanced SQM in the fig- 
ure signifies scenarios where both A and D are substantially 
elevated, causing the CFL SQM to be unable to maintain pres- 
sure equilibrium. This leads to the unstable existence of the 
SQS. 


IV. SUMMARY 


In this study, we investigated the GW echoes generated by 
SQSs within the framework of an equivparticle model. Dis- 
tinct from prior research that relied on the basic MIT bag 
model with predetermined sound velocities, our approach uti- 
lized EOSs for SQM in the equivparticle model enriched with 
density-dependent quark masses. This integration not only 
encapsulates confinement but also quark pairing effects. Sig- 
nificantly, our findings emphasized the crucial role of quark 
pairing effects to enable a photon sphere and thereby facilitate 
the production of a GW echo. As the value of A increases, the 
EOS for SQM becomes more stiff, allowing the mass-radius 
relation of SQS to intersect with the photon sphere line and, 
therefore, yield GW echoes. Additionally, for specific model 
parameters A and D, the SQS bearing the maximum mass 
showcases the lowest GW echo frequency. We also estimated 
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